\(\int (a+\frac {b}{x})^p (c+\frac {d}{x})^q \, dx\) [269]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (warning: unable to verify)
   Maple [F]
   Fricas [F]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 19, antiderivative size = 96 \[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=-\frac {b \left (a+\frac {b}{x}\right )^{1+p} \left (c+\frac {d}{x}\right )^q \left (\frac {b \left (c+\frac {d}{x}\right )}{b c-a d}\right )^{-q} \operatorname {AppellF1}\left (1+p,-q,2,2+p,-\frac {d \left (a+\frac {b}{x}\right )}{b c-a d},\frac {a+\frac {b}{x}}{a}\right )}{a^2 (1+p)} \]

[Out]

-b*(a+b/x)^(p+1)*(c+d/x)^q*AppellF1(p+1,2,-q,2+p,(a+b/x)/a,-d*(a+b/x)/(-a*d+b*c))/a^2/(p+1)/((b*(c+d/x)/(-a*d+
b*c))^q)

Rubi [A] (verified)

Time = 0.05 (sec) , antiderivative size = 96, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.158, Rules used = {382, 142, 141} \[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=-\frac {b \left (a+\frac {b}{x}\right )^{p+1} \left (c+\frac {d}{x}\right )^q \left (\frac {b \left (c+\frac {d}{x}\right )}{b c-a d}\right )^{-q} \operatorname {AppellF1}\left (p+1,-q,2,p+2,-\frac {d \left (a+\frac {b}{x}\right )}{b c-a d},\frac {a+\frac {b}{x}}{a}\right )}{a^2 (p+1)} \]

[In]

Int[(a + b/x)^p*(c + d/x)^q,x]

[Out]

-((b*(a + b/x)^(1 + p)*(c + d/x)^q*AppellF1[1 + p, -q, 2, 2 + p, -((d*(a + b/x))/(b*c - a*d)), (a + b/x)/a])/(
a^2*(1 + p)*((b*(c + d/x))/(b*c - a*d))^q))

Rule 141

Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^(p_), x_Symbol] :> Simp[(b*e - a*f
)^p*((a + b*x)^(m + 1)/(b^(p + 1)*(m + 1)*(b/(b*c - a*d))^n))*AppellF1[m + 1, -n, -p, m + 2, (-d)*((a + b*x)/(
b*c - a*d)), (-f)*((a + b*x)/(b*e - a*f))], x] /; FreeQ[{a, b, c, d, e, f, m, n}, x] &&  !IntegerQ[m] &&  !Int
egerQ[n] && IntegerQ[p] && GtQ[b/(b*c - a*d), 0] &&  !(GtQ[d/(d*a - c*b), 0] && SimplerQ[c + d*x, a + b*x])

Rule 142

Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^(p_), x_Symbol] :> Dist[(c + d*x)^
FracPart[n]/((b/(b*c - a*d))^IntPart[n]*(b*((c + d*x)/(b*c - a*d)))^FracPart[n]), Int[(a + b*x)^m*(b*(c/(b*c -
 a*d)) + b*d*(x/(b*c - a*d)))^n*(e + f*x)^p, x], x] /; FreeQ[{a, b, c, d, e, f, m, n}, x] &&  !IntegerQ[m] &&
 !IntegerQ[n] && IntegerQ[p] &&  !GtQ[b/(b*c - a*d), 0] &&  !SimplerQ[c + d*x, a + b*x]

Rule 382

Int[((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.), x_Symbol] :> -Subst[Int[(a + b/x^n)^p*((c +
 d/x^n)^q/x^2), x], x, 1/x] /; FreeQ[{a, b, c, d, p, q}, x] && NeQ[b*c - a*d, 0] && ILtQ[n, 0]

Rubi steps \begin{align*} \text {integral}& = -\text {Subst}\left (\int \frac {(a+b x)^p (c+d x)^q}{x^2} \, dx,x,\frac {1}{x}\right ) \\ & = -\left (\left (\left (c+\frac {d}{x}\right )^q \left (\frac {b \left (c+\frac {d}{x}\right )}{b c-a d}\right )^{-q}\right ) \text {Subst}\left (\int \frac {(a+b x)^p \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^q}{x^2} \, dx,x,\frac {1}{x}\right )\right ) \\ & = -\frac {b \left (a+\frac {b}{x}\right )^{1+p} \left (c+\frac {d}{x}\right )^q \left (\frac {b \left (c+\frac {d}{x}\right )}{b c-a d}\right )^{-q} F_1\left (1+p;-q,2;2+p;-\frac {d \left (a+\frac {b}{x}\right )}{b c-a d},\frac {a+\frac {b}{x}}{a}\right )}{a^2 (1+p)} \\ \end{align*}

Mathematica [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(206\) vs. \(2(96)=192\).

Time = 0.46 (sec) , antiderivative size = 206, normalized size of antiderivative = 2.15 \[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=\frac {b d (-2+p+q) \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q x \operatorname {AppellF1}\left (1-p-q,-p,-q,2-p-q,-\frac {a x}{b},-\frac {c x}{d}\right )}{(-1+p+q) \left (-b d (-2+p+q) \operatorname {AppellF1}\left (1-p-q,-p,-q,2-p-q,-\frac {a x}{b},-\frac {c x}{d}\right )+x \left (a d p \operatorname {AppellF1}\left (2-p-q,1-p,-q,3-p-q,-\frac {a x}{b},-\frac {c x}{d}\right )+b c q \operatorname {AppellF1}\left (2-p-q,-p,1-q,3-p-q,-\frac {a x}{b},-\frac {c x}{d}\right )\right )\right )} \]

[In]

Integrate[(a + b/x)^p*(c + d/x)^q,x]

[Out]

(b*d*(-2 + p + q)*(a + b/x)^p*(c + d/x)^q*x*AppellF1[1 - p - q, -p, -q, 2 - p - q, -((a*x)/b), -((c*x)/d)])/((
-1 + p + q)*(-(b*d*(-2 + p + q)*AppellF1[1 - p - q, -p, -q, 2 - p - q, -((a*x)/b), -((c*x)/d)]) + x*(a*d*p*App
ellF1[2 - p - q, 1 - p, -q, 3 - p - q, -((a*x)/b), -((c*x)/d)] + b*c*q*AppellF1[2 - p - q, -p, 1 - q, 3 - p -
q, -((a*x)/b), -((c*x)/d)])))

Maple [F]

\[\int \left (a +\frac {b}{x}\right )^{p} \left (c +\frac {d}{x}\right )^{q}d x\]

[In]

int((a+b/x)^p*(c+d/x)^q,x)

[Out]

int((a+b/x)^p*(c+d/x)^q,x)

Fricas [F]

\[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=\int { {\left (a + \frac {b}{x}\right )}^{p} {\left (c + \frac {d}{x}\right )}^{q} \,d x } \]

[In]

integrate((a+b/x)^p*(c+d/x)^q,x, algorithm="fricas")

[Out]

integral(((a*x + b)/x)^p*((c*x + d)/x)^q, x)

Sympy [F]

\[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=\int \left (a + \frac {b}{x}\right )^{p} \left (c + \frac {d}{x}\right )^{q}\, dx \]

[In]

integrate((a+b/x)**p*(c+d/x)**q,x)

[Out]

Integral((a + b/x)**p*(c + d/x)**q, x)

Maxima [F]

\[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=\int { {\left (a + \frac {b}{x}\right )}^{p} {\left (c + \frac {d}{x}\right )}^{q} \,d x } \]

[In]

integrate((a+b/x)^p*(c+d/x)^q,x, algorithm="maxima")

[Out]

integrate((a + b/x)^p*(c + d/x)^q, x)

Giac [F]

\[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=\int { {\left (a + \frac {b}{x}\right )}^{p} {\left (c + \frac {d}{x}\right )}^{q} \,d x } \]

[In]

integrate((a+b/x)^p*(c+d/x)^q,x, algorithm="giac")

[Out]

integrate((a + b/x)^p*(c + d/x)^q, x)

Mupad [F(-1)]

Timed out. \[ \int \left (a+\frac {b}{x}\right )^p \left (c+\frac {d}{x}\right )^q \, dx=\int {\left (a+\frac {b}{x}\right )}^p\,{\left (c+\frac {d}{x}\right )}^q \,d x \]

[In]

int((a + b/x)^p*(c + d/x)^q,x)

[Out]

int((a + b/x)^p*(c + d/x)^q, x)